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Abstract 

The three simultaneous algebraic equations, = 1, [C,VT\ = 0, [C,H] = 0, which determine 
the C operator for a non-Hermitian "PT-symmetric Hamiltonian H, are shown to have a nonunique 
solution. Specifically, theC operator for the Hamiltonian H = ip^ + i/x^g^+ieg^ is determined per- 
turbatively to first order in e and it is demonstrated that the C operator contains an infinite number 
of arbitrary parameters. For each different C operator, the corresponding equivalent isospectral 
Dirac-Hermitian Hamiltonian h is calculated. 
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I. INTRODUCTION 



A Hamiltonian H defines a physical theory of quantum mechanics if (i) H has a real energy 
spectrum, and (ii) the time-evolution operator U = e~*^* is unitary so that probability is 
conserved. These two features of the theory are guaranteed if H is Dirac Hermitian. (The 
Hamiltonian H is Dirac Hermitian ii H = H\ where the symbol f represents the combined 
operations of complex conjugation and matrix transposition.) However, it is not necessary 
for H to be Dirac Hermitian for the spectrum to be real and for time evolution to be unitary. 
For example, the Hamiltonians belonging to the class 

H = p^ + q\iqy (e>0) (l) 

possess real eigenvalues [l|, [2| and generate unitary time evolution [s], 0] . These Hamiltonians 
therefore define physically acceptable quantum theories. 

The Hamiltonians in ([1]) are VT symmetric, that is, invariant under combined spatial 
reflection V and time reversal T. The underlying reason that the Hamiltonians H in are 
physically acceptable is that they are selfadjoint, not with respect to the Dirac adjoint f, 
but rather with respect to CVT conjugation, where C is a linear operator that represents a 
hidden reflection symmetry of H. The CVT inner product defines a positive definite Hilbert 
space norm. Not every PT-symmetric Hamiltonian has an entirely real spectrum but if the 
spectrum is entirely real, a linear PT-symmetric operator C exists that obeys the following 
three algebraic equations [3|: 

C = 1, (2) 

[C,PT] = 0, (3) 

[C, H] = 0. (4) 

When such a C operator exists, we say that the VT symmetry of H is unbroken. Constructing 
the C operator is the key step in showing that time evolution for a non-Hermitian VT- 
symmetric Hamiltonian ([1]) is unitary. 

There has been much research activity during the past decade on non-Hermitian VT- 
symmetric Hamiltonians and the C operator for some nontrivial quantum-mechanical 
models has been calculated perturbatively 0, 0, S, 0] • The first perturbative calculation of 
C was performed in Ref. Q for the Hamiltonian 

H = \p' + + ieq\ (5) 

Here, /i is a mass parameter, e is a small coupling constant that is used as a perturbation 
parameter, and p and q are the usual quantum-mechanical operators satisfying the Heisen- 
berg commutation relation = i. It was shown in Ref. [6| that the C operator has a 
simple and natural form as the parity operator V multiplied by an exponential of a linear 
Dirac Hermitian operator Q: 

C = e^V, (6) 
where Q = ■ For the Hamiltonian Q is a series in odd powers of e: 

Q=Qie+ Qse' + Qse' + ■ • • ■ (7) 
To first order in e the result for Q was found to be 
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Note that in the unperturbed hmit e — in which the Hamiltonian becomes Hermitian and 
parity symmetry is restored, the operator Q vanishes and thus C — ^ P. This suggests that 
the C operator may be interpreted as the complex extension of the parity operator V. 

Mostafazadeh showed that the Q operator can be used to construct a similarity trans- 
form that maps the non-Dirac-Hermitian Hamiltonian H io a. spectrally equivalent Dirac- 
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h = e2/2i/e-2/2. (9) 

This similarity transformation was originally used by Scholtz et al. to convert Hermitian 
Hamiltonians to non-Hermitian Hamiltonians 

In Sec. HI] of this paper we reinvestigate the Hamiltonian i7 in ([5]). We show that the 
perturbative solution in ([8]) for the Q operator is in fact not unique and that the general 
solution to the algebraic system ([2]) - (j4]) contains an infinite number of arbitrary contin- 
uous parameters. Then, in Sec. IIIIl we show that for each of these C operators there is a 
corresponding spectrally equivalent Hermitian Hamiltonian h. 



II. COMPLETE FIRST-ORDER CALCULATION OF THE C OPERATOR 

The purpose of this paper is to show that the "PT-symmetric Hamiltonian in (j5]) has an 
infinite class of associated C operators. To determine the C operator for a given Hamiltonian 
H we must solve the three equations ([2]) - (jl]). The first two of these equations are kinematic 
constraints on C; they are obeyed by the C operator for any PT-symmetric Hamiltonian. If 
we seek a solution for C in the form ([6]), we find that ([2]) and ([3]) imply that q) is an odd 
function of the p operator and an even function of the q operator. The third equation (jl]) 
is dynamical because it makes explicit reference to the Hamiltonian. Our specific objective 
here is to solve (jlj) perturbatively to first order in e. 

We substitute ([5]) and ([6]) into (jlj) and to first order in e obtain the commutation relation 
satisfied by Qi: 

[Q^,Ho] = 2^q^ (10) 

where Hq = ^p"^ + ^f^^Q^ is the unperturbed Hamiltonian. The simplest solution to (fTOl) that 
is odd in p and even in q and is a polynomial in p and q is given by the solution in ([H]). Note 
that ([8]) is not the unique function of p and q that satisfies the commutation relation (fTOj) 
because we can add to Qi any function of Hq. However, Hq is an even function of p and Qi 
is required to be an odd function of p. Thus, one is tempted to conclude [wrongly]) that the 
only solution to ffTOl) is that given in ([8]). 

To construct additional solutions to ffTOj) we employ the operator notation described in 



detail in Ref. [12|. We introduce the set of totally symmetric operators where Tm,n is 
an equally weighted average over all possible orderings of m factors of p and n factors of q. 
For example, 

2o,o = 1) 
Tifi = p, 

= lipq + qp), 
Ti,2 = \{pqq + qpq + qqp), 

^2,2 = \ (ppqq + qqpp + pqqp + qppq + qpqp + pqpq) , (11) 
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and so on. Note that Tm,n is the quantum-mechanical generahzation of the classical product 
p^q^. We can express the polynomial solution Qi in ([8]) to the commutation relation (fTOj) 
as a linear combination of two such totally symmetric operators: 

The advantage of the totally symmetric operators T^^n is that it is especially easy to eval- 
uate commutators and anticommutators. For example, as shown in Ref. [l2|, the operators 
Tm,n obey extremely simple commutation and anticommutation relations: 

\Qi Tfji,n\ ^^-^m— 1,71) 

{(l)Tm^n\ = 2Tm,n+l- (13) 

Also, by combining and iterating the results in (fT3|) . we can establish additional useful 
commutation relations for Tm,n- For example, 

[g^^m,n] = 2i'mTrr,-l,n- (14) 

The totally symmetric operators Tm.n can be re-expressed in Weyl-ordered form (l2| : 

= ^ E (Dp'^y-' = ^ E (I) - = 0, 1, 2, 3, . . . ). (15) 

The proof that the totally symmetric form of Tm,n in ffTTj) equals the binomial-summation 
Weyl-ordered forms above requires the repeated use of the Heisenberg algebraic property 
that [q,p] = i and follows by induction. The reason for introducing the Weyl-ordered form 
of Tm,n is that it allows us to extend the totally symmetric operators Tm,n to negative values 
of n by using the first of these formulas or to negative values of m by using the second of 
these formulas. The commutation and anticommutation relations in f[T51) and ([T^ remain 
valid if m is negative or if n is negative. 

We will now show that (IT^ is just one of many solutions to flTUl) that are even in q 
and odd in p. Note that (fTOj) is an inhomogeneous linear equation for Qi and that (fT2|l 
is a particular solution to this equation. To find other solutions, we need only solve the 
associated homogeneous linear equation 

, homogeneous ) Ho] = 0. (16) 
Let us look for a solution to this equation of the form 

CO 

Qi , homogeneous 

fc=0 

which by construction is odd in p and even in q. Substituting (fT71) into (fT6|) and using the 
commutation relations in (fT4|) . we obtain a two-term recursion relation for the coefficients 

afc+i = --^Afc (A; = 0, 1, 2, 3, ■■■)• (18) 
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The solution to this recursion relation is 

= ;7^!^(z}y(-^')' (/c = 0, 1,2,3, ■■■), (19) 

where a is an arbitrary constant. Thus, a one-parameter family of solutions to the homoge- 
neous equation (fT6|) is 

Ql, homogeneous = — , / ly "^3-2fc,2fc- (20) 

We have included the factor of in the denominator so that eQi^ homogeneous is dimensionless. 

Combining the inhomogeneous solution in f|T2l) and the homogeneous solutions in fl20|) 
gives a general one-parameter class of solutions to the commutator condition ffTOj) : 

Qi = -^T3,o - ^Ti,2 + ^ g l^,r(-|) ( V)' T3_2M.. (21) 

The coefficients in this series are indicated schematically by the diagonal line of a's in Fig. [1] 
that intersects the m axis at (3, 0) and the n axis at (0, 3). On this figure we use circles to 
indicate the two terms in ( l2Tl) that arise from the particular polynomial solution (fT2l) . 

One might be concerned that because the operator p in the series in ( 1211) is raised to 
arbitrarily high negative powers it might be difficult to interpret the series as an operator. 
However, we will now show how to sum this operator series to obtain a well defined and 
meaningful result. First, we observe that there is an extremely simple way to represent the 
totally symmetric operator T_„ „ for all integers n: 

T-n,n = ]:{q-\ + . (22) 



2 \ pj 2 \p 

This nontrivial formula can be verified by using the Heisenberg algebra = i. One may 
regard this formula as the companion to the result that 

Tn,n = Hahn„ (Ti^i) , (23) 

where Hahn„(x) is a Hahn polynomial of degree n [l3| . 

In order to use the formula ( !2T|) . we replace the kth term in the series in (pT!) with a 
triple anticommutator of T_2k,2k'- 

T3.2k,2k = l{{{T.2k,2k,p},p},p}. (24) 

We then recognize that the sum is just a pair of binomial expansions of the general form 



oo 



k=0 ^ ' 

with a = |, one for each of the two operators and ^q. Thus, the one-parameter family 
of solutions in (!2T!) simplifies to 

Qi = -Ar3,o-ATi,2 + T^<;<;<;(i + /xVg-) .v).v).v 



3//^ ' /i2 ' 16/i' 





(26) 
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Note that at the classical level, where q and p commute, the p's in the anticommutators 
combine with the expression to the power | to give Hq^^, which is even in p. However, in 
(|26l) we can see that the symmetry requirement that Q be odd in p and even in q is satisfied! 
Furthermore, there are just enough powers of p in the triple anticommutators to cancel any 
small-p singularities. 

One might now be tempted to argue as follows: In the Heisenberg algebra a commutator 
behaves as a derivative. Thus, the commutator equation ffTOj) is analogous to a first-order 
linear ordinary differential equation. Since the solution to such an equation contains one 
arbitrary parameter, the solution in fl26l) . which contains the arbitrary parameter a, should 
be the complete solution to ffTOl) . However, this argument is wrong and we will now show 
that there are in fact an infinite number of additional one-parameter families of solutions to 

m- 

We seek new one-parameter families of solutions to the associated homogeneous com- 
mutation relation flTBl) . with each family of solutions labeled by the integer P = 
0, ±1, ±2, ±3, • ■ ■ . For each value of P the solution has the form 

oo 

2l'!homogcncous = '^r^^2P+3-2fc,2fc, (27) 

k=0 

which is odd in p and even in q, as is required. We substitute (127|) into (fT6!) and use the 
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FIG. 1: Schematic representation of the coefficients in the operator series solution for Qi in (j3ip . 
The two circles represent the terms in the inhomogeneous polynomial solution in (jl2p . The a's rep- 
resent the coefficients in the operator series ()20p , which is the homogeneous solution corresponding 
to P = 0. The coefficients of the homogeneous solutions in (|3Up for P = —1, —2, and —3 are 
indicated by 6's, c's, and d's, and the coefficients in the homogeneous solutions for P = 1, 2, and 
3 are indicated by x's, y's, and z's. 
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commutation relations f|T^ to obtain a recursion relation for the coefficients a\f^: 

ai% = - ^~^~K 'ar (fc = 0, 1, 2, 3, ■ ■ ■ ). (28) 
The solution to this recursion relation is 

"^'^ j^ Lr(lpl{! '-^''' (fc = 0.1.2,3,...), (29) 

where a^^^ are arbitrary constants. Thus, for each integer P we obtain the following one- 
parameter family of solutions to the homogeneous equation flTB]) : 

o(-) - ,2.^r (30) 

^1, homogeneous " 7TF7 P SY / ^ 2P+3-2fc,2fc • (.OUJ 

^ fc=0 ^--^ 2j 

The factor of /i^"^^ in the denominator ensures that eQi^homoeeneous dimensionless. Note 
that if we set P = 0, we recover the homogeneous solution in ll20|) . In Fig. [T]the coefficients 
for the P = —1, P = —2, and P = —3 solutions are indicated by fe's, c's, and d's and the 
coefficients for the P = 1, P = 2, and P = 3 solutions are indicated by x's, y's, and z's. 
We can now combine all these solutions to produce a very general solution for Q\ : 

^_ 4 2 >^ a(^) ^r(A;-F-|) ^ 

2l - -^^3,0 - -Jl,2 + 1^J^,1_,^ k\T{-P-l) ^"^ ^ ^2P+3-2.,2.. (31) 

This solution generalizes and replaces that in (!2T!) and is one of the principal results in this 
paper. We emphasize that a*^^'* are all arbitrary. 

Finally, we observe that it is possible to present the infinite sums over operators in (I5U]1 
and fl31l) more compactly by performing the sum over k. To do so, we generalize the result 
in flMl) as a 2P + 3-fold anticommutator for the case when P > — 1, 



^2P+3-2fe,2fc — q2P+z{{ - ■ ■ {T-2k,2k,P}, ■ ■ •P},P}2P+3 times • 

(32) 

This result allows us to apply the identity in to perform the sum in fl5Ul) . We find that 

^ ^xP+3/2 



^1, homogeneous 2"^^+^^!^+^ 1 1 ' ' ' I V P pj 




,pr ^■■■P? ^P 



2P+3 times 




22P+4^p+4 I I • • • M - ^ -^-^ ) ^P}^---P}^P} ■ (33) 



2P+3 times 

One may simplify this expression even further by exploiting the analog of the Baker- 
Hausdorff formula for multiple anticommutators 14| : 

e^Be^ = B + {A,B} + ^{A, {A, B}} + ^{A, {A, {A, B}}} + ■■■. (34) 

Thus, the 2P -|- 3-fold anticommutator in fl32|) . for example, may be rewritten in a more 
compact form as a multiple derivative: 



(^2P+3 

{{•••{7'_2fc,2fc,P},---P},P}2P+3timcs = ^^2P+3 '^^^ ^-2fc.2fc 



(35) 

/3=0 

When P < — 1, we can use commutators instead of anticommutators to make the expressions 
more compact, but we do not present the results here because they are repetitive. 
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III. SPECTRALLY EQUIVALENT HAMILTONIANS 



For the C operators associated with Qi in fl3Tl) we must now calculate the Hermitian 
Hamiltonian h that is equivalent (in the sense that it is isospectral) to the Hamiltonian in 
H in To do so, we evaluate the similarity transformation in which to leading order 
in e amounts to evaluating the commutator 



h 



Ho + \ie^ 



[q\ Qi] 



(36) 



Note that when we use the first- order form of the Q operator, we obtain h to second order 
in e. (The third-order contribution to Q gives h to fourth order in e.) 
In order to evaluate the commutator in fl36l) we use the relation 



(37) 



[q^, Tra,n\ = -jim{m - l)(m - 2)T„_3,„ + 3imTm-i,n+2, 

which is derived from the commutation relations in ( fT3l) and dH]). Furthermore, since (l36ll 
depends on Qi linearly, we can evaluate each of the contributions to Qi independently and 
add together the results at the end of the calculation. 

First, we evaluate the commutator in fl36|) for the inhomogeneous part (the first two 
terms) in Qi and obtain 



h = Ho + ^^{-l + m,2 + 3^2To,4) 



(38) 



This result was already obtained in Ref. [15|], where it was observed that to this order in per- 
turbation theory the equivalent Hermitian Hamiltonian h represents a quartic anharmonic 
oscillator having a position-dependent mass. However, this interpretation holds only for the 
case in which we take the coefficients a*-'^-* to vanish for all P. 

For each nonzero value of a^^^ the contribution to the commutator in ( l36l) for Q inhomogeneous 
in (I30D is 



homogeneous 



le a 



E 



^ ^ I) / .2\k r 3 m ] 

(~A* ) [Q 1 J 2P+3-2fc,2fcJ 



fc=0 



2/i 



P+4 



E 

fc=0 



2\k 



k\ 



[p + i-k)T{k-p + \) 3fcr [k-p-D ^ 

-l2P-2k,2k TTiTl 5^-'2P+4-2/fc,2A: 



r(-p-l) 



^^r (-P - f ) 



(39) 



where we have made use of ( !37l) . 

We have shown in this paper that the C operator contains an infinite number of arbitrary 
parameters. It remains an open question as to whether there is an additional physical or 
mathematical condition that would determine the C operator uniquely; that is, whether there 
is an advantageous or a "best" choice for the arbitrary parameters a^^^ in fl3Tl) . For example, 
one might try to choose a^^^ such that the equivalent Hermitian Hamiltonian h has the form 

+ V{q). [The non-Dirac-Hermitian "PT-symmetric Hamiltonian H = — is spectrally 
equivalent to the Dirac- Hermitian Hamiltonian h = p"^ + 4:gq^ — 2y/ghq [l^, 17, [H, 19, 20|.] 
However, there appears to be no choice of the parameters a^^^ that achieves such a simple 
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form. Thus, we conclude by posing the question, is there a useful fourth constraint that one 
can use to supplement the three conditions in ([2]), ([3]), and (jl]), that gives rise to the "best 
form" for the C operator? 

We conjecture that the answer to this question in the context of quantum field theory 
is that the huge parametric freedom in the C operator can be used to impose the condition 
of locality. A C operator for a g(p^ quantum field theory has been calculated perturbatively 
to first order in 51 0]. The resulting Q operator was found to decay exponentially at large 
spatial distances because it contains an associated Bessel function. It may be possible to 
exploit the parametric freedom in C to replace the Bessel function by a spatial delta function 
so that the C operator becomes local. 
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